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Abstract 



o 

A multiple exp-function method to exact multiple wave solutions of nonlinear partial dif- 
ferential equations is proposed. The method is oriented towards ease of use and capability of 
computer algebra systems, and provides a direct and systematical solution procedure which gen- 
eralizes Hirota's perturbation scheme. With help of Maple, an application of the approach to 
the 3 + 1 dimensional potential- Yu-Toda-Sasa-Fukuyama equation yields exact explicit 1-wave 
and 2-wave and 3-wave solutions, which include 1-soliton, 2-soliton and 3-soliton type solutions. 
Two cases with specific values of the involved parameters are plotted for each of 2-wave and 
3-wave solutions. 

PACS codes: 02.30.Gp, 02.30.1k, 02.30.Jr 
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O ■ 1 Introduction 

7— I 

o 

Exact solutions to nonlinear partial differential equations help us understand the physical phe- 
nomena they describe in nature. Many solution methods have been proposed, which contain the 
tanh-function method [H [21 [3] , the sech-function method [HOE], the homogeneous balance method 
[3[5] the extended tanh-function method j2J-[H], the sine-cosine method [HI 03], the tanh-coth 
method [2] and the exp-function method |15} [T6] . The crucial idea of these methods is to search 
for rational solutions to variable coefficient ordinary differential equations transformed from given 
nonlinear partial differential equations. Following this observation, a unified approach to exact 
solutions to nonlinear equations has been proposed, revealing relations between solvable ordinary 
differential equations and nonlinear partial differential equations recently [17] . Solitary waves, pe- 
riodic waves and kink waves modeling various nonlinear motions have been presented for many 
nonlinear dispersive and dissipative equations, indeed. 
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However, those existing methods are only concerned about travelling wave solutions to nonlinear 
equations. It is known that there are multiple wave solutions to nonlinear equations, for instance, 
multi-soliton solutions to many physically significant equations including the KdV equation and the 
Toda lattice equation [18J and multiple periodic wave solutions to Hirota bilinear equations [l~9] l2"0]. 
Therefore, it naturally comes that there should be a similar direct approach for constructing multiple 
wave solutions to nonlinear equations. We would, in this paper, like to give an answer by formulating 
a solution algorithm for computing multiple wave solutions to nonlinear equations. The approach 
will be illustrated step by step while applying to an example, providing a general feature of solving 
nonlinear equations by adopting linear ones. 

The application example we will present is the 3 + 1 dimensional so-called potential- Yu-Toda- 
Sasa-Fukuyama equation (for short, the potential- YTSF equation): 

- 4u xt + u xxxz + 4u x u xz + 2u xx u z + 3u yy = 0. (1.1) 

This equation is a potential-type counterpart of a 3 + 1 dimensional nonlinear equation 

[-4v t + $(v)v z ] x + 3v yy = 0, $ = d 2 + Av + 2v x d-\ (1.2) 

introduced by Yu, Toda, Sasa and Fukuyama in |21j . while making a 3+1 dimensional generalization 
from the 2 + 1 dimensional Calogero-Bogoyavlenkii-Schiff equation (see, say, |22j and references 
therein): 

- Av t + §(y)v z = 0, $ = d 2 + Av + 2v x d~ 1 , (1.3) 

as did for the KP equation from the KdV equation. Taking v = u x transforms the equation (|1.2p 
into the potential- YTSF equation (jl.ip [23]. We also remark that the equation (jl.ip itself becomes 
the potential KP equation if z = x, and reduces to the potential KdV equation while further taking 
Uy = 0. Therefore, various applications of the KP and KdV equations show great potential for 
applications of (jl.ip in the physical sciences. 

Obviously, the potential- YTSF equation (jl.ip has the solutions independent of two variables: 

u = f{z,t), u = f(x) + g(t), u = cx + f(z), u = cy + f(z), u = cy + f(t), (1.4) 

and a particular variable separated solution: 

u = (cy + d)x + yf(z,t) + g(z,t), (1.5) 

where c, d are arbitrary constants and /, g, h are arbitrary functions in the indicated variables; and 
a known solution u = u(x, y, z, t) will lead to a new one: 

v = u(x,y,z,t) + cy + f(t), (1.6) 

where c is an arbitrary constant and / is an arbitrary function in t. Moreover, a Backlund trans- 
formation of the type v = 2(\n(j)) x + u was constructed by Yan in [23] and a class of other variable 
separated solutions was constructed in |24]-|27], It is worth noting that variable separated solutions 
exist ubiquitously for 2 + 1 dimensional integrable equations (see, say, |28j). We will formulate 
a multiple exp-function solution method and present a few broad classes of exact wave solutions, 
including 1-soliton, 2-soliton and 3-soliton type solutions, to the potential- YTSF equation (jl.ip . 
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In particular, our multiple exp-function method will yield two different classes of two-wave and 
three-wave solutions to the potential- YTSF equation, and every class contains diverse soliton type 
solutions, both analytic and singular. 

The paper is organized as follows. In Section[2l a direct formulation for generating multiple wave 
solutions to nonlinear equations is established, by searching for rational solutions in new variables 
defining individual waves. In Section [3j an application is made to construct multiple wave solutions 
to the 3 + 1 dimensional potential- YTSF equation. We conclude the paper in the final section, along 
with a discussion on polynomial solutions. 

2 A multiple exp-function method 

Let us formulate our solution procedure by focusing on a scalar 1 + 1 dimensional partial differential 
equation 

P(x,t,u x ,u t ,~') = 0, (2- 1 ) 

which is assumed to be of differential polynomial type like the KdV equation. The solution method 
will also work for systems of nonlinear equations and high-dimensional ones. 

Step 1 - Defining Solvable Differential Equations: 

We introduce a sequence of new variables rji = rji(x, t), 1 < i < n, by solvable partial differential 
equations, for instance, the following linear ones: 

m,x = kVi, Vi,t = -^iVi, 1 < i < n , ( 2 - 2 ) 

where k{, 1 < i < n, are the angular wave numbers and Wj, 1 < i < n, are the wave frequencies. 
This is often a starting point for constructing exact solutions to nonlinear equations, since no way 
can help solve nonlinear equations directly. Solving such linear equations leads to the exponential 
function solutions: 

r\i = CiS\ £j = kiX — ojit, 1 < i < n, (2-3) 

where q, 1 < i < n, are any constants, positive or negative. The arbitrariness of the constants 
(H, 1 < i < n, brings more choices for solutions than we used to |29j . Each of the functions r]i, 

1 < i < n, describes a single wave and a multiple wave solution will be a combination of all those 
single waves. We emphasize that the linear differential relations in (|2.2p are extremely helpful while 
transforming differential equations to algebraic equations and carrying out related computations by 
computer algebra systems. The explicit solutions (|2.3p offer reasons why the approach is called the 
multiple exp-function method. The idea of using linear differential conditions could also be applied 
for other occasions, in which there might be diverse solutions |30j . Both the differential relations 
and the solution formulas are important in understanding and applying the approach. 

The basic idea of using solvable differential equations was also successfully used to solve the 

2 + 1 dimensional KdV-Burgers equation through a second-order ordinary differential equation 
ar]" + brj' + erf + dr] = (a,b,c,d = const.) in [31], and the Kolmogorov-Petrovskii-Piskunov 
equation through a first-order ordinary differential equation r/ = 1 ± rj 2 in [9] . It has been broadly 
adopted in the tanh-function type methods [lOj [HI Q3] , the Jacobi elliptic function method [32| [33] . 
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the mapping method [34J|35], the F-expansion type methods [361 [371 [38] an d the G'/G-expansion 
method [39]. 



Step 2 - Transforming Nonlinear PDEs: 

Let us proceed to consider rational solutions in the new variables rji, 1 < i < n: 

p(rj 7] T] ) n M n N 

Q[Vl,V2,- ■■ ,Vn) r , s = UJ= o r,s=li,j=0 

where Pknj and qkUj are ah constants to be determined from the original equation (|2.ip . All Laurent 
polynomial and polynomial functions are only special examples of rational functions, and so, we 
can similarly have a multiple tanh-coth method for getting multiple wave solutions to nonlinear 
equations. 

By using the differential relations in (j2.2j) , it is straightforward to express all partial derivatives 
of u with x and t in terms of r/j, 1 < i < n. For example, we can have 

n n n n 

9 Y PViVht 9*^.* ~ q Y U iP^Vi +PY "ilViVi 

i=l i=l i=l i=l /r> r \ 

U t = 5 = 2 ' ( 2 - 5 ) 

and 

n n n n 

q PvrVi,x -pY QviVi,* q ^PviVi -pY k ^viVi 



i=l i=l i=l i=l 



q2 q2 



(2.6) 



where p Vi and q m are partial derivatives of p and q with respect to r/j. This way, we can see that 
all partial derivatives, not only ut and u x , will still be rational functions in the new variables r]i, 
1 < % < n. Substituting those new expressions of partial derivatives into the original equation (|2.ip 
generates a rational function equation in the new variables rji, 1 < i < n: 

Q(x,t,m,V2, ■■■ ,Vn) = o. (2.7) 

This is called the transformed equation of the original equation (|2.ip . The step here makes it 
possible to compute solutions to differential equations directly by computer algebra systems. 

Step 3 - Solving Algebraic Systems: 

Now we let the numerator of the resulting rational function Q(x, t, rji, rj2, ■ ■ ■ ,rj n ) to be zero. 
This yields a system of algebraic equations on all variables ki,uji,pki ; ij,qki,ij] an d solve this system 
to determine two polynomials p and q and the wave exponents £j, 1 < i < n. All computation can 
be done systematically by computer algebra systems such as Maple. We point out that the resulting 
algebraic systems may be complicated and so a computer program really helps. Now, the multiple 
wave solution u is computed and given by 

u(x,t) g( Cl gfeia;— k>i* ; . . . ; c n e knX ~ u)nt ) ' (^"^) 

Since we begin with the exponential function solutions to the initial linear equations, we call the 
above method a multiple exp-function method. If we choose some other linear equations, we can, 
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for instance, have a multiple sine-cosine method to get multiple periodic wave solutions to nonlinear 
equations. Clearly, our multiple exp-function method in the case of n = 1 becomes the so-called 
exp- function method proposed by He and Wu in [15j . 

The solution procedure described above provides a direct and systematical solution procedure 
for generating multiple wave solutions and it allows us to carry out the involved computation 
conveniently by powerful computer algebra systems such as Maple, Mathematica, MuPAD and 
Matlab. It also presents a generalization of Hirota's perturbation scheme to construct multi-soliton 
solutions [18] . We will analyze three cases of polynomials p and q for the 3 + 1 dimensional potential- 
YTSF equation (jl.ip . to construct its multiple wave solutions. 



3 One-wave, two-wave and three-wave solutions to the potential- 
YTSF equation 

Let us apply our multiple exp-function method to the 3 + 1 dimensional potential- YTSF equation 
(|1.1|) . We will discuss three cases of two polynomial functions p and q to generate one- wave, two- 
wave and three- wave solutions as follows. 

Case 1 - One-wave solutions: 

We require the linear conditions: 

m,x = hm, vi, y = km., m,z = m^, m,t = -vim., (3.1) 

where ki,k,mi,ui are constants. Then try a pair of two polynomials of degree one: 

p(m) = a o + a ivi, q(vi) = bo + hm, (3-2) 

where ao, at, &o> b\ are constants to be determined. By the multiple exp-function method and using 
the differential relations in (|3.ip , we can have the following solution to the resulting algebraic system 
with Maple: 

6i(2A;i& + a ) 1. 2 3 k 2 fo , 

ax = Fq , = --h mx - — , (3.3) 



and all other constants are arbitrary. Since we can have an exponential function solution to 

rji = e k i x + l iy+ m i z ~ UJ i t ; (3-4) 

the corresponding 1-wave solutions read 

, . p a + a^+hy+m.z-^t 

U = U{X,y, Z.tj = - = : : : 7, (3.5) 

v ' y ' ' ' g 5 Q _|_ l) ie kix+hy+miz-ujit ' v I 

where a\ and uj\ are defined by (|3.3|) and all the other involved constants are arbitrary. This is in 
agreement with the selection for the 1-soliton solution in [30] and contains all exact solutions in 
[41] . Note that the wave frequency depends on all angular wave numbers in the 1-wave solutions 
above, but we will see that it is not the case in the 2-wave and 3-wave solutions below. 

Case 2 - Two-wave solutions: 
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Similarly, we require the linear conditions: 

Vi,x = kVii Vi,y = kVU m,z = mVh Vi,t = -WW, 1 < % < 2, (3.6) 
where ki,li,mi,uJi, 1 < i < 2, are constants, and thus, the solutions 771 and 772 can be defined by 

m = ^kix+hy+miz-wit ^ 1 < j < 2 . (3.7) 

where ci and C2 are arbitrary constants. 

Let us try a particular pair of two polynomials of degree two: 

p(vi,m) = 2[hm + k 2m + aMh + £2)771772], 

(3.8) 

Q{V1>V2) = 1 + ??i + ??2 + ai 2 rnV2, 

where ai2 is a constant to be determined. By the multiple exp-function method and using the 
differential relations in ()3.6p . we can have two solutions to the resulting algebraic system with 
Maple: 

3 1 

uii = —-^h — -ki 2 m,i, 1 < i < 2, (3-9) 

and 

when li = ki, 1 < ? < 2; and 

1 3/- 2 

c^ = -^ 3 -^p 1 <7<2. (3.11) 

and 

(fcifc 2 2 - ki 2 k 2 + fci? 2 ~ hk2) (kifo 2 ~ ki 2 k 2 - hh + hk 2 ) 
° 12 ~ {hk 2 2 + /ci 2 £ 2 + ^1^2 - hk 2 ) {k x k 2 2 + fci 2 ^ - k x l 2 + ^£2) ' 
when rrii = ki, 1 < i < 2. 

Then, the two corresponding 2-wave solutions are determined by 

, ,n P{m,m,m) ^[hm +k 2 r] 2 + a 12 (k x +k 2 )rn 772] , > 

u = u x, y, z, t) = — = — ■ ■ , 3.13) 

q{riuV2,m) 1 + m + m + a i2ViV2 



where 771 and 772 are defined by (|3.7p . either with the frequencies uj\ and uj 2 being given by (|3.9p 
and 012, by (|3.10p when Zj = ki, 1 < i < 2; or with the frequencies ui\ and uj 2 being given by 
(|3.1ip and 012, by (|3.12p when mi = ki, 1 < i < 2. All the unspecified involved constants in the 
solutions are arbitrary. There is a different selection of frequencies in [12] but it does not lead to 
exact non-constant solutions. Two specific solutions of the above 2-wave solutions are plotted in 
the figures IBTTI and I3T2"1 In each figure, the first plot is three dimensional, and the other plots exploit 
the x-, y- and z-curves or the contour plots with z = at different times. 

Case 3 - Three-wave solutions: 

Again similarly, we require the linear conditions: 

r]i,x = ki-rji, rji >y = kr]i, r] i)Z = rrnrji, 77^ = -w^j, 1 < i < 3, (3.14) 

where ki,li,rrii,uji, 1 < i < 3, are constants, and thus, the solutions 771,772 and 773 can be defined by 

m = c-ekix+iiv+miz-uit j 1 < • < 3 ( 3 _ 15 ) 
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Figure 3.1: 1st 2-wave solution with k% = 1, k 2 = —2, m\ = 1, m 2 = 5, c\ = 1, c 2 = 2. 



where ci,C2 and C3 are arbitrary constants. 

Let us now try the following particular pair of two polynomials of degree three: 
P(f7i,%,?73) = 2[/cir/i + k 2 r] 2 + k 3 r] 3 + ai 2 (A:i + k 2 )r]im + ai 3 (&i + ^3)^1% 

+a 23 (£; 2 + k 3 )r] 2 ri 3 + ai 2 ai 3 a 23 (/ci + & 2 + A; 3 )?7i r? 2 r/ 3 ], (3.16) 

Qim^m^m) = ^ + m + m + m + avmxm + 01^1% + 023%% + ai 2 ai 3 a 23 r/ir? 2 ?73, 

where ai 2 , ai 3 and a 23 are constants to determined. By the multiple exp-function method and using 
the differential relations in (|3.14p . we can have two solutions to the resulting algebraic system with 
Maple: 



and 



when li = ki, 1 < i < 3; and 



3 1 

LOi = --h - -ki 2 m,i, 1 < i < 3, (3.17) 

aij= ti+ k S ' i -^'- 3 ' (3 - i8) 

1 3/ 2 

^ = -^^ 3 -^, 1 < i < 3, (3.19) 



1 < < 3, (3.20) 



and 22 22 

^k^kj ~\~ k% kj ~\~ kjl j Jj^kj^ ^k^kj ~\~ k^ kj kjlj ~\~ l^kj^ 
when rrii = ki, 1 < i < 3. 

Then, the two corresponding 3-wave solutions are given by 

u = u(x, y, z, t) = — -, (3.21) 

where p and q are defined by (|3.16|) and 771 , ry 2 and r/ 3 are defined by (|3.15|) , either with the frequencies 
ui,u>2 and a; 3 being given by (|3.17f) and a\ 2 , ai 3 and a 23 , by (|3.18p when li = ki, 1 < i < 3; or 
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Figure 3.2: 2nd 2-wave solution with k\ = 1, k2 = 3, l\ = 2,l<i = 1, c\ = 1, c% = 2. 



with the frequencies uii,U2 and U3 being given by (|3.19p and 012, 013 and £123, by (|3,20p when 
rrii = ki, 1 < % < 3. All the unspecified involved constants in the solutions are arbitrary. Two 
specific solutions of those 3-wave solutions are plotted in the figures 13.31 and 13.41 In each figure, 
the first plot is three dimensional, and the other plots exploit the x-curves with y = 1 and different 
z- values at different times or the contour plots with z = at different times. 

We emphasize that through the proposed multiple exp-function algorithm, two kinds of 2-wave 
solutions and 3-wave solutions are easily obtained for the potential- YTSF equation (jl.lj) . If for 
2-wave and 3-wave solutions, we take the general wave frequencies like (j3.3jl . where mi, ki,l\ have 
no relation, we will meet contradictions in the resulting algebraic systems. On the other hand, if the 
involved constants in (|3.5p satisfy &o^i < and some of the constants Cj, 1 < i < n, in (|3.13p and 
(|3,2ip are negative, the corresponding exact solutions become singular. Moreover, for the second 
case (i.e., rrii = ki), even if the constants Cj are positive in (|3.13|) and (|3.2ip . the constants aij can 
be negative, and thus, the solutions (|3. 13|) and (|3.2ip can be singular. Taking special constants in 
our 1-wave, 2-wave and 3-wave solutions and considering equal angular wave numbers = rrii = ki 
yields all special soliton solutions to the potential-YTSF equation (jl.ip . presented by Wazwaz in 
03]. 

4 Concluding remarks 

A direct and systematical solution procedure for constructing multiple wave solutions to nonlinear 
partial differential equations is proposed. The presented method is oriented towards ease of use 
and capability of computer algebra systems, allowing us to carry out the involved computation 



8 




-250 -200 — 150 — 100 —SO O -250 -200 — 1 50 — 1 OO — 50 O 



Figure 3.3: 1st 3-wave solution with k% = 0.8, k 2 = 1.6, k 3 = — 0.6, /i = — 2, Z2 = 3, Z3 = — 1.5, ci = 
0.9, c 2 = 0.8, c 3 = 1.2. 

conveniently through powerful computer algebra systems. It is the use of computer algebra systems 
that in each case of 2-wave and 3-wave solutions, we are able to present two classes of concrete 
exact explicit solutions to the 3 + 1 dimensional PYTSF equation, only in form of u = f(t, x + y,z) 
or u = f(t,x + z,y) (but not in a general form including u = f(t,x,y + z)). The key point of 
our approach is to search for rational solutions in a set of new variables defining individual waves. 
An application of our method yields specific 1-wave, 2-wave and 3-wave solutions to the 3 + 1 
dimensional PYTSF equation. The method can also be easily applied to other nonlinear evolution 
and wave equations in mathematical physics. 

It is direct to check that the 3 + 1 dimensional potential- YTSF equation has the following 
class of polynomial solutions: 

ui = ui(x, y, z, t) = a\ + a 2 x + a 3 y + a 4 z + a 5 t + a G xy + a^yz + a 8 yt + agzt + a w yzt, (4.1) 

where a^, 1 < % < 10, are arbitrary constants. These are all polynomial solutions among a class 
of polynomial functions with deg(«i,x) = deg(iti,y) = deg(iti,z) = deg(ni,t) = 1. On the other 
hand, there are other two solutions: 

u 2 = u 2 {x, y, z, t) = ai + a 2 x + a 3 y + a^z + a 5 t + f(z, t), (4.2) 

and 

n 3 = u 3 (x, y, z, t) = a\ + a 2 x + a 3 z + a 4 t + g(x) + h(t), (4.3) 

where aj, 1 < i < 5, are arbitrary constants and f,g,h are arbitrary functions in the indicated 
variables. Taking /, g, h as polynomials engenders other polynomial solutions to the potential- 
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Figure 3.4: 2nd 3-wave solution with k% = 0.8, k<i = 1.6, = — 0.6, ii = —2,1% = 3,^3 = — 1.5, ci = 
0.9, c 2 = 0.8, c 3 = 1.2. 

YTSF equation (|1.1|) . which can be of high degree. But the third one reduces to solutions to the 
2 + 1 dimensional potential Calogero-Bogoyavlenkii-Schiff equation, independent of the variable y. 

It is our guess that higher-wave solutions to the 3 + 1 dimensional potential- YTSF equation 
(|1.1|) could be presented in a parallel manner. But the required computation is pretty complicated, 
even in the case of 4-wave solutions. We hope that they could be presented and verified by some 
analytic way. Any general form of 2-waves and 3-waves, which does not involve any relation among 
the angular wave numbers ki,li,rrii, will be more interesting and important. 
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